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Abstract 

We study multidimensional gravitational models with scalar curvature nonlinearities of the type and 
i?*. It is assumed that the corresponding higher dimensional spacetime manifolds undergo a spontaneous 
compactification to manifolds with warped product structure. Special attention is paid to the stability of 
the extra-dimensional factor spaces. It is shown that for certain parameter regions the systems allow for 
a freezing stabilization of these spaces. In particular, we find for the _R~^— model that configurations with 
stabilized extra dimensions do not provide a late-time acceleration (they are AdS), whereas the solution 
branch which allows for accelerated expansion (the dS branch) is incompatible with stabilized factor spaces. 
In the case of the i?'* —model, we obtain that the stability region in parameter space depends on the total 
dimension D = dim(M) of the higher dimensional spacetime M. For D > 8 the stability region consists of 
a single (absolutely stable) sector which is shielded from a conformal singularity (and an antigravity sector 
beyond it) by a potential barrier of infinite height and width. This sector is smoothly connected with the 
stability region of a curvature-linear model. For _D < 8 an additional (metastable) sector exists which is 
separated from the conformal singularity by a potential barrier of finite height and width so that systems in 
this sector are prone to collapse into the conformal singularity. This second sector is not smoothly connected 
with the first (absolutely stable) one. Several limiting cases and the possibility for inflation are discussed 
for the J?*— model. 



PACS numbers: 04.50.-Hh, 11.25.Mj, 98.80.Jk 

1 Introduction 

Distance measurements of type la supernovas (SNe,Ia) [1] as well as cosmic microvifave background (CMB) 
anisotropy measurements [2] performed during the last years give strong evidence for the existence of dark 
energy — a smooth energy density with negative pressure which causes an accelerated expansion of the Universe 
at present time. This late time acceleration stage should have started approximately 5 billion years ago and 
represents a second acceleration epoch after inflation which lasted for 10""^^ seconds immediately after Big Bang 
and ended 13.7 billion years ago. 

The challenge to theoretical cosmology consists in finding a natural explanation of inflation and late time 
acceleration (dark energy) within the framework of string theory/M-theory or loop quantum gravity. Scenarios 
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which address one or both of these issues arc, e.g., string (pre-bigbang) cosmology [3], a large number of brane 
world scenarios (brane inflation [4] , the ekpyrotic [5] and born-again Universe [6] scenario, late-time acceleration 
via S-branes [7]), as well as the string theory/M-theory scenarios with flux compactifications [8, 9, 10, 11, 12] 
and the recent setups of loop quantum cosmology [13] . 

A different starting point for explaining the late time acceleration was taken in Refs. [14, 15], where it 
was shown that purely gravitational modifications of the (curvature linear) Einstein-Hilbert action by including 
curvature nonlinear terms of the type R~^, could induce a positive effective cosmological constant, and 
with it an accelerated expansion. The corresponding phenomenological studies were performed in four spacetime 
dimensions by transforming the curvature-nonlinear theory into an equivalent curvature-linear theory with the 
nonlinearity degrees of freedom carried by an additional dynamical scalar field. (The used technique had been 
developed in earlier work since the 1980s [16, 17, 18, 19].) A natural question which arises with regard to this 
phenomenological approach is of whether it can naturally follow as low-energy limit from some M-theory setup 
or — looking from down-side up — of whether it has a physically viable phenomenological extension to higher 
dimensions. 

In the present paper we take the latter (phenomenological) point of view and study higher dimensional 
extensions of purely gravitational non-Einsteinian models with scalar curvature nonlinearities of the types R^^ 
and R'^. Special emphasis will be laid on finding parameter regions (regions in moduli space) which ensure the 
existence of at least one minimum of the effective potential for the volume moduli of the internal spaces and which 
by this way allow for their stabilization. The latter fact is of special importance because the extra-dimensional 
space components should be static or nearly static at least from the time of primordial nucleosynthesis (otherwise 
the fundamental physical constants would vary, sec e.g. [20, 21], and observational bounds on the variation 
of the fine-structure constant could be violated [23]). During the last few years, problems of volume moduli 
stabilization have been studied, e.g., for models with large extra dimensions (Arkani-Hamed-Dimopoulos- 
Dvali (ADD) setups [24]) [25], as well as (more recently) for M-theory scenarios with flux compactifications 
[9, 10, 11, 12, 27, 28] and for brane-gas systems [29]. 

Here wo will mainly follow our earlier work on this subject [23, 30, 31, 32, 33, 34], where the stabilization 
of extra dimensions was studied for {Dq + Z)')— dimensional bulk spacetimes with a product topology. The 
corresponding product manifolds are constructed from Einstein spaces Mi with scale (warp) factors which de- 
pend only on the coordinates of the external Dq— dimensional spacetime Mq (the ansatz resembles a zero mode 
approximation in Kaluza-Klein formalism). As a consequence, the excitations of the scale factors (conformal ex- 
citations/excitations of the volume moduli) have the form of massive scalar fields (gravitational excitons/radions) 
living in the external spacetime. Stabilized volume moduli will correspond to positive eigenvalues of the mass 
matrix of these scalar fields, unstable configurations to tachyonic excitations. 

The present work can be understood as a direct continuation of our earlier investigation on volume moduli sta- 
bilization in i?^— models of purely geometrical type [33] as well as with magnetic (solitonic, Preund-Rubin-type) 
form fields living in the extra dimensions (flux fleld compactifications) [34] . Its key results can be summarized 
as follows: 

• A straight forward extension of a four-dimensional purely geometrical _R^^— model to higher dimensions 
with subsequent dimensional reduction cannot simultaneously provide a late time acceleration and a 
stabilization of the extra dimensions. A late time acceleration is only possible for a solution branch which 
has a positive definite maximum of the effective potential, i.e. a positive definite effective cosmological 
constant, and not a negative definite minimum as it would be required for a stabilization of the internal 
factor space components. This means that other, more sophisticated, extension scenarios would be needed 
to reach both goals simultaneously. 

• In contrast to the .R^— models of Refs. [33, 34], the .R^— model shows a rich substructure of the stability 
region in parameter (moduli) space which crucially depends on the total dimension D = Dq + D' of the 
bulk spacetime. There exists one stability sector which is present for all dirnensions D > Dq + 2 and 
which smoothly tends to the stability sector of an _R— linear model when the _R— nonlinearity is switched 
off. This sector is shielded from the (probably unphysical) antigravity sector of the theory by a potential 
barrier of infinite height and width, and hence, it will be absolutely stable with regard to transitions of the 
system into the antigravity sector. Apart from this sector of absolute stability, there exists a second sector 
for total dimensions D < 8 which is separated from a conformal singularity and an antigravity sector 
beyond it by a potential barrier of finite height and width. Systems in this sector will only be mctastable 
and prone to collapse into the conformal singularity and the antigravity sector. The metastable sector is 
separated from the stability sector of the i?— linear model by an essential singularity in R and the effective 
potential Ueff (when R and Uefj are considered as functions over a parameter subspace). 

First discussions of this subject date back to Refs. [22]. 
2 See also [26]. 
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• In the analyzed purely geometrical _R— nonlinear models, a stabilization of the internal factor spaces is 
necessarily connected with an AdS structure of the external spacetime. In order to uplift the AdS to an 
dS with positive effective cosmological constant and the possibility for late time acceleration, additional 
matter fields should be included into the model (e.g., flux fields as in the ii^— model of Ref. [34]). 

The paper is structured as follows. In section 2, we present a brief technical outline of the transformation 
from a non-Einsteinian purely gravitational model with general scalar curvature nonlinearity of the type f{R) 
to an equivalent curvature- linear model with additional nonlinearity carrying scalar field. Afterwards, we derive 
in Sec. 3 criteria which ensure the existence of at least one minimum for the effective potential of the internal 
space scale factors (volume moduli). These criteria are then used in sections 4 and 5 to obtain the regions 
in parameter (moduli) space which allow for a freezing stabilization of the scale factors in models with scalar 
curvature nonlinearities of type R^^ and R"^. The main results are summarized in the concluding Sec. 6. 



2 General setup 

We consider a £> = (4 + D')— dimensional nonlinear pure gravitational theory with action functional 



= 2^/ 



x^fiR) , (1) 

where f{R\ is an arbitrary smooth function with mass dimension 0{m?') (m has the unit of mass) of a scalar 
curvature R = R[g\ constructed from the D— dimensional metric gab {o-,b = 1, • • • , D). D' is the number of extra 
dimensions and k|, denotes the D— dimensional gravitational constant which is connected with the fundamental 
mass scale M*(4+£3/) and the surface area Sd-i = 27r(-^~^^/^/r[(L' — l)/2] of a unit sphere in Z) — 1 dimensions 
by the relation [35, 36, 37] 

Kl = 2SD-i/M'^^,^'^,y (2) 
The equation of motion for the theory (1) reads (see, e.g., Refs. [16, 17, 18]) 

f'Rab -\f9ab- ^a^bf + Qab^f = (3) 

and has as trace 

(D - 1)0/' = |/ - f'R . (4) 

We use the notations Vo and □ for the cqyariant derivative and the Laplacian with respect to the metric Qab, 
as well as the abbreviations /' = df/dR, Rab = Rabig]- 

Before we endow the metric of the pure gravity theory (1) with explicit structure, we recall that this 
_R— nonlinear theory is equivalent to a theory which is linear in another scalar curvature R but which contains 
an additional self- interacting scalar field. According to standard techniques [16, 17, 18], the corresponding 
i?— linear theory has the action functional: 



(5) 



where 



and where the self-interaction potential [/((/>) of the scalar field cf) is given by 

u{4>) = lifr''^^''~'^[Rf'-f] , (7) 

= ie-^n i?(</.)e^'^ - / (i?(0))] , B := ^ ^ (8) 

The metrics Qab, Sab and the scalar curvatures R, R of the two theories (1) and (5) are conformally connected 
by the relations 

gab = n^gab=[nR)f^''~^^gab (9) 
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and 



D-2' 



(10) 



via the scalar field (j) = ln[f'{R)]/A. This scalar field (f) carries the nonlinearity degrees of freedom in R of the 
original theory, and for brevity we call it the nonlinearity field. 

As next, we assume that the D-dimensional bulk space-time M undergoes a spontaneous compactification 
to a warped product manifold 

M = Mo X Ml X . . . X Mn (11) 

with metric 



(12) 



i=l 



The coordinates on the (.Do = rfo + 1)— dimensional manifold Mq (usually interpreted as our observable [Dq 
4)— dimensional Universe) are denoted by x and the corresponding metric by 



(13) 



For simplicity, we choose the internal factor manifolds Mj as dj— dimensional Einstein spaces with metrics 
5^'^ = 9mini{yi)dy^' (S) dy^', so that the relations 



Rn 



rrii, rii 



l,...,di 



and 



R 



X^di — Ri 



(14) 



(15) 



hold. The specific metric ansatz (12) leads to a scalar curvature R which depends only on the coordinates x of 
the external space: R\g\ = R{x). Correspondingly, also the nonlinearity field (j) depends on x only: (j) — 4'{x)- 

Passing from the i?— nonlinear theory (1) to the equivalent iJ— linear theory (5) the metric (12) undergoes 
the conformal transformation g ^ g [see relation (9)] 



n 



2r(a:)„(i) 



with 



g(0) 



(3' := f3' 



A 



D-2' 



(16) 
(17) 



3 Freezing stabilization 

The main subject of our subsequent considerations will be the stabilization of the internal space components. 

A strong argument in favor of stabilized or almost stabilized internal space scale factors /3* (x) , at the present 
evolution stage of the Universe, is given by the intimate relation between variations of these scale factors and 
those of the fine-structure constant a [23]. The strong restrictions on a— variations in the currently observable 
part of the Universe [39] imply a correspondingly strong restriction on these scale factor variations [23]. For 
this reason, we will concentrate below on the derivation of criteria which will ensure a freezing stabilization of 
the scale factors. Extending earlier discussions of rnodels with R^ scalar curvature nonlinearities [33, 34] we will 
investigate here models of the nonlinearity types R~^ and R"^. 

In Ref. [32] it was shown that for models with a warped product structure (12) of the bulk spacetime M 
and a minimally coupled scalar field living on this spacetime, the stabilization of the internal space components 
requires a simultaneous freezing of the scalar field. Here we expect a similar situation with simultaneous freezing 
stabilization of the scale factors /3*(.t) and the nonlinearity field (i)(x). According to (17), this will also imply a 
stabilization of the scale factors /?' (x) of the original nonlinear model. 

In general, the model will allow for several stable scale factor configurations (minima in the landscape over 
the space of volume moduli). We choose one of them, denote the corresponding scale factors as (}q, and work 
further on with the deviations 

(x) = f3' (x) - f5i (18) 

^For a discussion of possible decompactification scenarios we refer to the recent work [38]. 
Although the toy model ansatz (1) is highly oversimplified and far from a realistic model, we can roughly think of the chosen 
minimum, e.g. , as that one which we expect to correspond to the current evolution stage of our observable Universe. 
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as the dynamical fields. After dimensional reduction of the action functional (5) we pass from the intermediate 
Brans-Dicke frame to the Einstein frame via a conformal transformation 



-2/(Do-2) 



3(0) 



(19) 



with respect to the scale factor deviations /3' {x) [33, 34, 36] . As result we arrive at the following action 
1 



S = 



2/^ 7-1 

^0 



(20) 



Mo 



which contains the scale factor offsets through the total internal space volume 

VD'=ViXVo = f[[ d^'yJ]: 

J Mi 



15^1 X I I e 



in the definition of the effective gravitational constant nj-,^ of the dimensionally reduced theory 

An 



^(00=4) = i^d/Vd' = Sn/Mi 



S 



D-l 



*(4+D')" 



(21) 



(22) 



Obviously, at the present evolution stage of the Universe, the internal space components should have a total 
volume which would yield a four- dimensional mass scale of order of the Planck mass M(4) = Mpi. The tensor 
components of the midisuperspace metric (target space metric on M^) G,j = 1, . . . ,n), its inverse metric 
G*-' and the effective potential are given as [40, 41] 



1 _.. 5»J 1 



The effective potential has the explicit form 



no-2 



i=l 



-2/3* 



(23) 



(24) 



where we abbreviated 

:=i?iexp(-2/3^). (25) 

For completeness, we note that the original metric gab of the ^—nonlinear model and the final Einstein frame 
metric g^i, of the dimensionally reduced model are connected by the relation 



9ub = {f) °- 



e2r^(.) 



(26) 



which up to the nonlinearity induced conformal factor (j')~^/(^~^) takes (for scale factors depending only on 
the time coordinate) a similar form like in the recently analyzed cosmological S'— brane models of Refs. [7, 42] . 

A freezing stabilization of the internal spaces will be achieved if the effective potential has at least one 
minimum with respect to the fields l3^{x). Assuming, without loss of generality, that one of the minima is 
located at /?* = /3o = 0, we get the extremum condition: 



eff 



dpi 



= 



Ri 



0=0 



Do 



-J2Rj + 2U{cI>) 



(27) 



^For simplicity, wc work here with stabihzed scale factors /3q which we assume as homogeneous and constant. In general, one 
can split the scale factors f3^{x), e.g., into a coherent scale factor background Pf){x) and non-coherent scale factor fluctuations 
/3*(x-) = 0^{x) — Pq{x) over this background [31]. 

''An alternative stabilization mechanism can consist, e.g., in the recently proposed dynamical stabilization in the vicinity of 
enhanced symmetry points [43]. In our present discussion we will not analyze such scenarios. 



5 



From its structure (a constant on the l.h.s. and a dynamical function of 4>{x) on the r.h.s) it follows that a 
stabilization of the internal space scale factors can only occur when the nonlinearity field 4>(x) is stabilized as 
well. In our freezing scenario this will require a minimum with respect to (f): 



= 



dU, 



eff 



d(i) 



= 0. 



(28) 



We arrived at a stabilization problem, some of whose general aspects have been analyzed already in Refs. 
[30, 32, 33, 34]. For brevity wo only summarize the corresponding essentials as they will be needed for more 

detailed discussions in the next sections. 

1. Eq. (27) implies that the scalar curvatures Ri and with them the compactification scales e^'o [see relation 
(25)] of the internal space components are finely tuned 



Ri 
di 



i,j = l,...,n. 



(29) 



2. The masses of the normal mode excitations of the internal space scale factors (gravitational excitons/radions) 
and of the nonlinearity field ^ near the minimum position are given as [32] : 



nil 



(i2 [/(</,) 



> 0. 



(30) 
(31) 



3. The value of the efii'ective potential at the minimum plays the role of an effective 4D cosmological constant 
of the external (our) spacetime Mq: 



eff 



/3»=0, 



(32) 



4. Relation (32) implies 



signAe// = sign[/((/)o) = signiJj 



(33) 



Together with condition (30) this shows that in a pure geometrical model stable configurations can only 
exist for internal spaces with negative curvature Ri < {i = 1, . . . ,n). Additionally, the effective cosmo- 
logical constant Ag// as well as the minimum of the potential U{(j)) should be negative too: 



Ae//<0, U{cPo)<0. 



(34) 



5. Eqs. (21), (22), (25) and (29) - (32) yield the following scaling behavior of the minimum related model 
parameters under a change of one of the offset scale factors /3o(i) — * (^o(2) ■~ ^^W(^o{i) ^ ^^^^ 
induced, e.g., by a change of the minimum value C/(0o(i)) — ^ C^(0o(2))- 



g/3S(2) 



2 



m 



fe(i) 



R 



*(i) 



A, 



m 



fe(2) 



R 



A, 



e//(2) 



''(Oo=4),(l) 



D'(2) 



K 



(no=4),(2) 



V, 



D'{1) 



1 D'/2 



U{(l>0{2)). 



(35) 
(36) 

(37) 



''Negative constant curvature spaces Mi are compact if they have a quotient structure: M, = H'^* /Ti, where H"^* and Fj axe 
hyperbolic spaces and their discrete isometry group, respectively. 
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6. For a system which is almost stabilized at a freezing point 
constant background and small deviations: 



Ao, the function f{R) can be split into a 



f{R) « ci{R - Ro) + /(^o) = ciR + C2 



(38) 



where ci ,f'{Ro) = exp{A(f)o), Ro = Ri^Pa), and — C2/(2ci) plays the role of a cosmological constant. 
In the case of homogeneous and isotropic spacetime manifolds, linear purely geometrical theories with 
constant A— term necessarijy imply_an anti-deSitter/deSitter geometry so that the manifolds are Einstein 
spaces. Substitution of f{R) ciR + C2 into Eq. (3) proves this fact directly 



Rab 



C2. 



£>- 2ci 



R 



D C2 



(39) 



Pliigging the potential U{(t>) from Eq. (8) into the minimum conditions (28), (31) yields with the help of 
d4,R = A/7/" the conditions 



dU 



A 



if 



-D/{D-2) 



h : 



4>o 



2{D - 2) 
Df-2Rf', ^ /^(0o)=O, 

1 , .,.-2/{D-2) 1 



(/': 



<t>0 



2{D-l) 

where the last inequality can be reshaped into the suitable form 

f"dRh\^, = f" [{D - 2)/' - 2Rf"]^^ > . 

Furthermore, we find from Eq. (40) 

D-2 



>0, 



C/(<^o) = 

so that (34) leads to the additional restriction 



2D 



(/')-^ R{^o) 



R{M < 



(40) 



(41) 



(42) 



(43) 



(44) 



at the extremum. Via the relation (10) the stabilization point curvatures of the i?— nonlinear and the i?— linear 
models are connected as 

Ro « (/')^/(^-^)i?o. (45) 

Thus, as the extra dimensional scale factors approach their stability position the bulk spacetime curvature 
asymptotically (dynamically) tends to a negative constant value (see Eq. (44)). Because the effective cosmo- 
logical constant is also negative (Ag// < 0), the homogeneous and isotropic external (Do — 4)— dimensional 
spacetime is asymptotically AdSuQ. Together with the compact hyperbohc internal spaces Mj = H'^^/Ti this 
results in a spontaneous compactification scenario 



AdS 



D 



AdSco X H'^'/Tix ...X H'^^ /r„ . 



(46) 



In the next sections we will analyze the conditions (29) - (41), (44) on their compatibility with particular 
scalar curvature nonlinearities f{R). 



4 The -model 

Recently it has been shown in Refs. [14] that cosmological models with a nonlinear scalar curvature term 
of the type R~^ can provide a possible explanation of the observed late-time acceleration of our Universe 
within a pure gravity setup. The equivalent linearized model contains an effective potential with a positive 
branch which can simulate a transient inflation-like behavior in the sense of an effective dark energy. The 
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corresponding considerations have been performed mainly in four dimensions. Here we extend these analyses 
to higher dimensional models — assuming that the scalar curvature nonlinearity is of the same form in all 
dimensions. We start from a nonlinear coupling of the type: 

fiR) = R-tx/R, /x>0. (47) 

In front of the _R^^— term, the minus sign is chosen, because otherwise the potential U{(f)) will have no extremum. 

With the help of definition (6), we express the scalar curvature R in terms of the nonlinearity field 4> and 
obtain two real- valued solution laranches 

E±=±VM(e^"^-l)"'^' , =^ '/'>0 (48) 
of the quadratic equation f'{R) = e^'^. The corresponding potentials 

U±{(l)) = ±^e-'^'t' Ve^"^ - 1 (49) 
have extrema for curvatures [see Eq. (40)] 



Ro,± = iVA* 



D 



D-2 



and take for these curvatures the values 



The stability defining second derivatives [Eq. (41)] at the extrema (50), 



-D/{D-2) 



D [DT2f 2D 



show that only the negative curvature branch _R_ yields a minimum with stable internal space components. 
The positive branch has a maximum with U+{4>o) > 0. According to (33) it can provide an effective dark energy 
contribution with A^ff > 0, but due to its tachyonic behavior with d'^U{(po) < it cannot give stably frozen 

internal dimensions. This means that the simplest extension of the four-dimensional purely geometrical R^^ 
setup of Refs. [14] to higher dimensions is incompatible with a freezing stabilization of the extra dimensions. 
A possible circumvention of this behavior could consist in the existence of different nonlinearity types fi{Ri) 
in different factor spaces Mi so that their dynamics can decouple one from the other. This could allow for a 
freezing of the scale factors of the internal spaces even in the case of a late-time acceleration with Ae/f > 0. 
Another circumvention could consist in a mechanism which prevents the dynamics of the internal spaces from 
causing strong variations of the fine-structure constant a. The question of whether one of these schemes could 
work within a physically realistic setup remains to be clarified. 

Finally we note that in the minimum of the effective potential Ueff{<t>, /?'), which is provided by the negative 
curvature branch R(4'), one finds excitation masses for the gravexcitons/radions and the nonlinearity field (see 
Eqs. (30), (51) and (52)) of order 

nil = ■ ■ ■ = rrin ^ ^ ^^^^ . (53) 

For the four-dimensional effective cosmological constant A^/f defined in (32) one obtains in accordance with 
Eq. (51) Ae// ~ 



5 The i?4-model 

In this section we analyze a model with curvature-quartic correction term of the type 

f{R)=R + jR^-2AD. (54) 

This setup contains no quadratic curvature terms and can be understood as a very rough approximative analogue 
of specific curvature corrected models of M-theory (see e.g. [19, 45, 46]). The investigation will be performed 

*A discussion of pro and contra of a higher dimensional origin of R^^ terms can be found in Kef. [15]. 
^The role of curvature-quartic corrections in M-theory inflation scenarios was recently discussed in Ref. [44]. 
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for an arbitrary number of dimensions, D. 

We start by deriving the explicit form of the potential U{4>). For this purpose we substitute f{R) from (54) 
into relation (6), 

/' = e^"^ = 1 + 47^^ (55) 

and resolve the latter equation for R: 

R = (47)- V3 (e^^ -if'^, -00 < (A < 00 . (56) 
The potential U {(p) is then found from Eq. (8) as 

(57) 



^(47)-^/^ {e^t'-lf' + 2Kn 



From its form with 

f/[?!>;Ar,>0,7>0]>0 V.;^ (58) 

we immediately conclude that the minimum condition U{<i>o) < cannot be satisfied in the sector (Ad > 
0, 7 > 0). In the remaining sectors, the potential will have a minimum if the extremum condition (40) and the 
minimum-ensuring inequality (42) will be fulfilled simultaneously. In the present case, these conditions read 

h[AD,i, R] ■■= b{D - 8)i?4 + {D- 2)R - 2DKd\ = (59) 

and 

f"df,h = 127i?2 [{D - 2) + A{D - 8)7i?-'] > 0, (60) 

respectively. From their structure it follows that the dimension D = % will constitute an exceptional class of 
models [due to cancellation of the highest order terms in (59), (60)]. We will analyze this class of models in 
subsection 5.2 below. 



5.1 Dimensions D 8 

Eq. (59), /i[Ai5,7, _R] = 0, is an algebraic equation in the variables {Ad,J,R) which defines a two-dimensional 
algebraic variety V C in the three-dimensional parameter space M. = M'^ 9 {Ad,^,R)- On this variety 
inequality (60) together with the restrictions (44) and (6) 



R<0, /' = 1 + 47^?'"' > 



(61) 



selects the parameter subset T C V of stably compactified internal space configurations. Choosing Ad and 7 as 
independent parameters, our main task will consist in obtaining the projection Q(AD,-y) '■— ''^'^ of the stability 

region T C V C onto the (A/), 7)— plane. (By tt we denote the projection itself.) Most of the_information 
will be derived by finding restrictions on (A£),7) from the conditions which ensure the reality of R as solution 
of the algebraic equation (59). 

In order to obtain the solutions R of equation (59) explicitly, we follow standard techniques (see, e.g., [47] 
and A) and consider first the associated cubic equation 



8DA 



D 



1{D- 



D-2 



(62) 



and its discriminant Q: 



Q = r^ + q\ 



q :-- 



8DAr 



3j(D - 



r :- 



D 



1 2 



-/{D - 8) 



(63) 



Depending on the sign of Q, the cubic equation (62) has one real solution for Q > or three real solutions for 
Q < 0, where in the case Q = at least two of these solutions coincide. Denoting (one of) the real solution(s) 
by Ml, the four roots of the quartic equation (59) can then be obtained according to (138), (139) as solutions of 
the two quadratic equations 

1 



R' 



±^R+- 



3q 







(64) 
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with 



-sign 



D-2 



(65) 



Physically sensible solutions will correspond to real roots of these equations. 

Following this general scheme of analysis, we start from the discriminant Q which we rewrite for later 
convenience as 



w{D) : 



(66) 



(£)-2)4 ■ 

In B it is show that the minimum ensuring inequality (60) implies 

0(Az,,7)>-l (67) 
so that Q is necessarily positive definite, Q > 0, and, hence, Eq. (62) has only one real-valued root [47] 



r^/^vi{z) > 0, 



)l/2 



1/3 



where 



vi{z) := l + {l + z) 



1/2 



1/3 



1-(1 + ^) 



1/2 



1/3 



(68) 



(69) 



It is now an easy task to explicitly analyze the pair of quadratic equations (64). They will only have real- valued 
roots, if at least one of the corresponding discriminants A± will be non-negative 



A± = -uiT'2e\Jul + 3q>0 . (70) 
(The subscripts ± in A± correspond to the signs in Eq. (64).) Because of Ui > this holds only for 

A_^ = -Ul+2^J^i4 + 3q (71) 
and under the additional reality-ensuring requirement 

ul + 3q>0. (72) 
Using the definitions (66), (69) of z and vi{z), the inequalities A_e > and (72) can be reshaped into the form 



Hi{z) := vf{z) + 6Az>0 
H2{z) := vf{z) + 27z>0, 



(73) 
(74) 



respectively. 

From the graphics of the functions Hi^^iz) depicted in Fig. 1, we read off that both inequalities (73), (74) 
are satisfied for ^ > —1 and no additional restrictions are set by them on the allowed region of the parameter z. 

Hence, we arrive at the result that the real- valued roots of the quartic equation (59) follow from the equations 
(64) by identifying in them ± = — e. The corresponding quadratic equation reads 



and has solutions (we distinguish them again by subscripts ±) 

1 



(75) 



R. 



e^/ul ± \ 2 J ul + 3q- Ui 



= ri/%±(0), 



(76) 
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Figure 1: Auxiliary functions Hx{z) and H^i^z), which for Hi^2{z) > ensure a non- negative discriminant 
A_g > of the quadratic equations (64). 



where 



{z) := \ l^^,± ^2^vf + 3z^/^-~?j 



(77) 



These roots are defined over the complete parameter region z{Ad,j) > —1 so that further restrictions on (A/j, 7) 
can only follow from the additional requirements (61). 

The first of these requirements, R < 0, should be fulfilled for a successful freezing stabilization of the extra- 
dimensional factor spaces. From the structure of (76), (77) we read off that contradicts this bound, 
whereas for the remaining solutions it partially narrows the allowed parameter region as follows 



R+,- 
R-,+ 

R_ _ 



< z, 

-1 < < 0, 
-1 < z. 



(78) 



The second inequality, /' 
restrictions 



1 + 47^?^ > 0, of (61) is analyzed in C and maps into the following parameter 



D < 
D > 



z < -w{D) = \w{D)\ , 
-w{D) < z. 



(79) 
(80) 



So far, wc have performed our analysis mainly in terms of the function 2; (A/), 7) and, hence, in terms of 
projections of the bounds (60) and /' > on the (A^i, 7)— plane. For completeness, we have to test whether all 
of the projected segments of V over the allowed region of the {Ad, 7)— plane fulfill the additional bound R <0, 
i.e. we should re-analyze (60) and /' > in terms of R. We start with (60). A simple case analysis gives 



7 > : 



7 < 



D > 
D< 

D > 
D< 



e = + 



e = + 



R^ > 
R^ < 

R^ > 
R^ < 



1 D-2 



1 D-2 



47£) - 8 
1 D-2 



47D-; 
1 D 



>0, 



47D-8 



(81) 
(82) 

(83) 
(84) 



^'^The restrictions (78) are only necessary conditions for the existence of solutions il < of the quartic equation, and provide 
only partial information about (60) and /' > 0. 
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and shows that configurations with Z? > 8, 7 < violate the bound R < 0, whereas (82) is weaker than R < 0. 
The remaining two inequalities (81), (84) can be reshaped with the help of (77) and 



1 D 



47D-8 



2-3/2^1/2 



(85) 



(from Eq. (63)) as 

respectively. From the graphics of the functions T__±(z) shown in Fig. 2 we read off that r_,+(z) > —2~-^/-^, 



0.5 




Figure 2: Rescaled scalar curvatures T-^±{z) := r ^^^R-^±. 

T^.-iz) < -2^1/2 foj. ^ > _i Hence, inequahty (81) is fulfilled only by R^^+, whereas (84) selects R. 
comparison of the inequalities (81) - (84) with the condition /' = 1 + 47_R'^ > and its implications 

7 > : -(47)-^ <R^, 7 < : i?^ < 1471-^ 



(87) 



shows that the latter relations (87) are compatible with them and add no additional restrictions. 

Combining the information about the inequalities (81) - (84) with (58) and the relations (78), (79), (80) 
we arrive at the following parameter regions for configurations with a possible freezing stabilization of extra- 
dimensional (internal) factor spaces: 



D > 8 : 
D<8: 



7 > 
7 < 
7 > 
7 < 



R-,+ , -w{D) < z <0, 
no stability , 

R+-, 0<z<\w{D)\, 
R--, -l<z<\w{D)\. 



(89) 
(90) 
(91) 



This result is schematically shown in Figs. 3,4. (The regions of formal extension to configurations with /' < 
are included in the graphics for reasons of completeness. They are briefly discussed in subsection 5.3, below.) 



5.2 The exceptional dimension D — 8 

In this particular case, the fl— nonlinear terms in the minimum-ensuring conditions (59), (60) cancel and the 
bounds on the solution can be read off immediately: 

R = ^Ad < 0, 7 > . (92) 

^'^For completeness we note that in the setup of the previous subsection 5.1 the formal limit — » 8 (for D assumed as non-discrete 
and real-valued) would correspond to the exceptional (singular) case 2 — » 0, Q — > -|-oo. 
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0.5 
>- 
-0.5 
-1 
-1.5 



/ 

z>-w(D) / 




f'<0 / f >o 1 
/ 




1— z = -w (D) 




^^^^ 0<z<-w(D)J 








-l<z<0 


0<z<-w(D) y"''^ 




z = -w (D)-/ 


z = -l A 


f'>0 /f'<0 



-0.4 



-0.2 




An 



0.2 



0.4 



Figure 3: Projection 0(Aj-,_^) of the stability region T C V C 7W of a subcritical model with D < 8 on the 
(A£),7)— plane (shaded areas with /' > 0). The two lines z(A£),7) = —w{D), given in Eq. (152), correspond 
to the conformal singularity f — where R —(47)^^/'^ maps into R —00. They separate parameter 
regions with U((f) —00) +00 (dark grey) from regions with U{(f> —00) —00 (light grey). For /' > 
the corresponding regions ensure the existence of an absolutely stable minimum (dark grey) and a metastable 
minimum (light grey, see the discussion in subsection 5.3 below). The (/' < 0)— regions have been included into 
the graphics for completeness. 



[R. < follows from relation (43) and condition (44).] The restriction /' = l + 47i?^ > sets the same additional 
bound on the allowed parameter region 



-1<4j(^^Ad^ <0 (93) 
as (80) in the case of D > 8. 

5.3 Stable and metastable configurations 

In the previous two subsections 5.1 and 5.2, it has been shown that models in D >8 dimensions posses a stability 
region 6(A£,.7) = ttT, which is located in the upper (A/j, 7)— half- plane. For models in £> < 8 dimensions this 
region is larger and extends additionally into the lower (A/j, 7)— plane (see the relations (88) - (91) and Figs. 
3, 4). Now we will analyze the system over these stability regions in more detail. 

We start with the asymptotic behavior of the potential U{(l>) in the limits (f> ±00. In the high curvature 
limit (j) +00 with dominant nonlinearity of the type /' = e"^'^ = 1-1- 47 ^ 1 we find from (57) 

U{ct> ^ +00) « |(47)-V3e(-B+4A/3)0 (94) 



and the sign of the coefficient in the exponent 

'-A-B = ^^ (95) 
3 3 D-2 ^ ' 

leads for fixed 7 to a qualitatively different behavior for dimensions D > 8 and D < 8: 

for £) > 8 , 

+00) ^ ^(47)-'/^ X <; 1 for£» = 8, (96) 

for D < 8 . 




^■^Hcrc and in the subsequent considerations the high curvature hmits are understood as formal hmits (in the mathematical 
sense) within the framework of our simplified toy model. It is clear that in a model with scalar curvature nonlinearity of the type 
/(-R) = R + '^j^^2'^kR'° gravitational self-interaction effects will dominate when 51^2'^*-^''"^ ~ ^- (^'^^ discussions of related 
subjects we refer to [48].) This means that a self-consistent treatment of the model would require techniques from (loop) quantum 
gravity or the high-energy sector of M-theory — what is out of the scope of the present paper. 
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z<-w (D) 
f ' <0 / 


f >o 1 

z=-w (D) 









-0.4 -0.2 0.2 0.4 



Ad 



Figure 4: Projection Q^Aon) of the stability region T C V C of a model with £) > 8 on the (A/j, 7)— plane 
(shaded areas with /' > 0). The line z{Ad,j) = —w{D), given in Eq. (152), corresponds to the conformal 
singularity /' = where R —(47)^^/''' maps into i? — 00. It separates the parameter region with 
U{(j) —00) +00 (dark grey) from the region with U{(f> —00) —00 (light grey). For /' > the 
corresponding region ensures the existence of an absolutely stable minimum (dark grey). The (/' < 0)— region 
is depicted for completeness. The lower half-plane 7 < corresponds to unstable configurations. 



The existence of a critical dimension (in our case -D = 8) is a rather general feature of gravitational theories 
with polynomial scalar curvature terms (see, e.g., Refs. [17, 49]). It can be easily demonstrated for a model 
with curvature nonlinearity of the type 

N 

f{R) = Y,akR' (97) 

fe=0 



for which the ansatz 



leads, similar like (7), to a potential 



N 



e^'^ = f = Y,kakF&-^ (98) 



fe=0 



N 



U{^) = \ (/')-^/^^-') Y.{k - l)a,RK (99) 

In the limit (f) — > +00 the curvature will behave like R « ce^'^ where h and c can be defined from the dominant 
term in (98): 

e^^ « Na^R''-' « iVa^vc^-'e^^-^)'*^. (100) 

Here the requirement /' > allows for the following sign combinations of the coeSicients ajv and the curvature 
asymptotics R{(j) 00): 

N = 21 : sign [a at] = sign [^(0 00)] 

TV = 2^ + 1: aN>0, sign [i?(0 -> 00)] = ±1. (101) 

The other combinations, N = 21 : sign [ajv] = —sign [R{(f) 00)], N = 2/-I-1 : oat < 0, sign [R{(l) 00)] = ±1, 
would necessarily correspond to the /' < sector, so that the complete consideration should be performed in 
terms of the extended conformal transformation technique of Ref. [18]. Such a consideration is out of the scope 

of the present paper and we restrict our attention to the cases (101). The coefficients h and c are then easily 

1 

derived as h = A/{N — 1) and c = sign(ajv) ]-^ctjv] ""^ . Plugging this into (99) one obtains 

U{(j) +00) « sign(ajv)^^^^ \NaN\~^ e'^^^e^^^^'^ (102) 
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and that the exponent 

D-2N , , , 

{D-2)iN-l) ^ (103) 

changes its sign at the critical dimension D = 2N: 



[/((/)-> +00) ^ sign (aAr)^——^ iTVaArp" X <^ 1 for D = 2N , (104) 



2N 




for D>2N, 



for D <2N. 



This critical dimension D = 2N is independent of the concrete coefficient ajv and is only defined by the degree 
dcgj^(/) of the scalar curvature polynomial /. From the asymptotics (104) we read off that in the high curvature 
limit (j) +00, within our oversimplified classical framework, the potential ?/((/)) of the considered toy-model 
shows asymptotical freedom for subcritical dimensions D < 2N, a stable behavior for ojv > 0, D > 2N and a 
catastrophic instability for qn < 0, D > 2N . We note that this general behavior suggests a way how to cure a 
pathological (catastrophic) behavior of polynomial nonlinear theories in a fixed dimension D > 2Ni: By 

including higher order corrections up to order iV2 > D /2 the theory gets shifted into the non-pathological sector 
with asymptotical freedom. More generally, one is even led to conjecture that the partially pathological behavior 
of models in supercritical dimensions could be an artifact of a polynomial truncation of an (presently unknown) 
underlying non-polynomial /(i?) structure at high curvatures — which probably will find its resolution in a 
strong coupling regime of M— theory or in loop quantum gravity. 




Figure 5: Typical form of a potential J7(0) with parameters (A^i, 7) from a subcritical {D < 8) region of absolute 
stability U{(j) —00) ^ +00, [/((/> — > -l-oo) — > -1-0 in the physical sector /' > 0. Specifically, it is set D — 
Ad = —1/4 for several values of 7 (in Planck units). 

As next step we consider the opposite limit cj) —00 which corresponds to /' — > 0. From (7) we see that 
the potential U{4>) of a model with general polynomial i?— nonlinearity (97) behaves in this limit like 

C/(^^_oo)«-ie-^^^/(i?c.), (105) 

where Rd is one of the real roots of the polynomial f'{R) — given in (98). In other words, the potential U{(j)) 
diverges like 

-00) -sign [fiRc^)] X 00 (106) 

for any dimension D and any value /(Rd) ^ 0. This means that, for energetic reasons, the system will 
be repelled from configurations with /' « by an infinitely high barrier U{4> —00) — > -|-oo in the case 
f{Rci) < 0, and it will be catastrophically attracted (experience a collapse) to /' w in the case f{Rci) > 0, 
U{(f) —00) — > —00. 

Let us make this general consideration explicit for the i?''— model (54). The polynomial /' = 1 + 47i?'^ has 
the single real-valued root Rc2 = —(47)"^^'^, which was used in (149) - (80) to map the inequality /' > into 
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Figure 6: Typical form of a potential U{(j)) of a subcritical {D < 8) metastable system with U{(t> ^ — oo) — > — cxa, 
[/((/> — > +00) ^ —0 in the physical sector /' > 0. It is set 13 = 6, ~ 1/4 for several values of 7 (in Planck 
units) . 



the (Ac, 7)— plane [result: the bounds ~w{D > 8) < z 

"3 



c2) 



(47)-i/3 



){D < 8)]. Plugging this root into / we get 



2A, 



From the condition /(i?c2) < for the existence of a repelling potential barrier with U{(j> —00) 
one finds the following inequalities (regardless of the existence of a minimum) 



(107) 
-00 limit 



7 > 
7 < 



^(47)i/^Az,>-l, 
|(47)i/=^Az,<-l. 



(108) 



(Below we will show that the case 7 > holds for the (/' > O)~sector, whereas 7 < will correspond to /' < 0.) 
Multiplication with w{D) leads to the equivalent conditions 



D>8: 
D <8: 



7 > 
7 < 
7 > 
7 < 



-w{D) < z, 

z < -w{D) , 

z < -w{D) = |w(i:>)| , 

\w{D)\ = -w{D) < z. 



(109) 
(110) 
(111) 

(112) 



We see that the asymptotical behavior of the potential U{4> — > —00) is defined by similar parameter regions in 
z like those which control the existence of a minimum under the condition /' > 0. In both cases the bound 
is connected with the critical value /' = which corresponds to Zc = —w{D). This is not a surprise because 
the two simultaneous conditions f'{Rc2) — 0, f{Rc2) — 0, which define the critical value Zc in the inequalities 
f{Rc2) < 0, f(Rc2) > 0, fulfiU the extremum condition d^U{(f>) = 0, i.e. h = Df - 2f'R = in (40) and hence 
the quartic equation (59). 

Analogously to inequalities (108) - (112) we get from the condition f{Rc2) > for a catastrophically 
attracting potential, U{(f> ^ —00) — > ^00, that such an asymptotics holds over the sectors 



D > 



D < 



7 > 
7 < 
7 > 
7 < 



z < -w{D) , 
-w{D) < z , 
\w{D)\ = -w{D) < z, 
z < ~w{D) = \w{D)\ . 



(113) 
(114) 
(115) 
(116) 



^Obviously, it holds trivially h = Df - 2f'R = for R = Rc2. 
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These sectors are eomplementary to those in (109) - (112). Here, we observe that 7 < for /' > and 7 > 
for /' < 0. This is confirmed by a comparison of inequaUties (109) - (116) with (79), (80) and the inequahties 
in 18. Obviously, the regions z < —w{D) for D > 8 and z > —'w{D) = \'w{D)\ for Z) < 8 correspond to a formal 
extension of the potential J7(^) into the (/' < 0)— sector. For completeness, these regions have been included 
into Figs. 3, 4. The typical U{<p — » —00) — » ±00 behavior of the potential is illustrated in Figs. 5,6. 

From Fig. 6 we see that the minimura is separated by a barrier of finite height and width from the singularity 
at ^ —00. Hence, we find that the i?**— models in the (/' > 0)— sector are absolutely stable for 7 > and 
metastable with tendency to collapse into the singularity /' w for 7 < 0. We further see from the conformal 
relation (10) between the scalar curvature R of the nonlinear model and the curvature R of the equivalent 
linear model that /' = 0, and hence z = \w{D)\, corresponds to a conformal singularity: The finite curvature 
value Rc2 = —(47)"^/^ of the nonlinear model is related to a curvature singularity R (/')~^^^^~^^-Rc2 in the 
associated linear model. 

A detailed study of various limiting cases is given in D. The corresponding main results can be summarized 
as follows. 

• The limit {Ajj —0,7 > 0) in the stable sector corresponds to a flat-space limit ^ ^ — which via 
(55) is associated with a freezing of the nonlinearity field (p: /' ^ 1 at (po 0. In the metastable sector 
nothing special happens in the limit (A^ — > 0, 7 < 0). 

• For Ad the limit 7 — > +0 corresponds to a freezing of the nonlinearity field at cpQ = and a 
smooth transition to a linear gravity model of Einstein-Hilbcrt type. In contrast, the limit 7 — 
of the metastable {D < 8)— system results in an infinitely deep minimum U{(j)o,J —0) —00 at 
(po = (IM) In \3D/{D - 8)1 and a curvature singularity ^-,-(7 ^ -0, £) < 8) w -\{D - 2)/['y{D - 8)]\^/^. 

• The limit z ^ — 1 (in the metastable sector) corresponds to coalescing minimum and maximum of the 
potential £/(</>) (with resulting inflection point at z = — 1). For z < —1 the potential U (0) has no minimum 
at all and the system is completely unstable. 



5.4 Inflation in the (7 > 0)— sector 

For simplicity we restrict our attention to most simple models with only one internal factor space Mi. The 
requirement for non-vanishing negative curvature of this space (in order to ensure a late-time stabilization of 
the corresponding dimensions) holds only for total dimensions > 6 (in the case of Dq =4). In terms of the 
normalized scale factor (radion) of this internal space Mi 



lch{D-2)^, 



Do 



(117) 



the effective potential (24) reads 



eff 



U{ct>) - -R,e^^^^ 



s :- 



(Do-2)(D-2)' 



£'0-2 

Si = -z S. 

di 



(118) 



In Figs. 7, 8 its generic form is illustrated by a model with {di = 2) extra dimensions and parameters 7 = 1/2, 
Ad = —1/4. A general feature of an eflFective potential (118) with a minimum and a barrier asymptotic 
U{4' —00) +00 is the necessary absence of a local maximum or a saddle point. This is easily seen from 
the extremum conditions 



t/(0) + -P°-J-^i fl,e^^i^' 



2di 



effl 



extr 



(119) 
(120) 



and their implications. Condition (120) yields again the quartic equation (59). For Z) > 8, 7 > this equation 
has only a single solution R-,+ in the (/' > 0)— sector so that the minimum is the only extremum. In the case 

^*Prom Fig. 2 and relation (88) one sees that the condition /' > cuts the T___|_(z)— branch at Zc2 = —w{D) (where —1 < Zc2 < 0) 
and allows only for the piece —w{D) < 2 < where T—^+{z) > T—^^{zc2) holds. The remaining segment —1 < 2 < —w{D) with 
T-^+{z) < T__+(2c2) and also the whole second solution T-^-(z) are completely located in the (/' < 0)— sector. 



17 



-10 



-4 -2 2 4 6 8 10 





Figure 7: Contour plot of the effective potential Ueff{(j),ip) given in Eq. (118). 



-D < 8, 7 > beside the negative curvature solution R+^- of the minimum there exists a positive maximum 
R+,+ > (see Eq. (76)) so that at the corresponding extremum with regard to (j) it holds U{(l)max) > in 
accordance with Eq. (43). Hence, the condition (119) cannot be fulfilled due to U{4)„iax) > 0, Ri < 0, Do = 
4, di > 2 =^ dipUeff > and there will be no other extremum of the effective potential Ueff{(t>,'i^) apart 
from the already studied minimum at (4> = (pQ, (p = 0). 

This means that in the considered oversimplified model inflation of purely topological type [50] as, e.g., 
recently demonstrated for SUGRA inspired setups (racetrack inflation starting at a saddle point of the effective 
potential) in Rcf. [51] is ruled out. In general, the too steep slopes of the exponential terms of the effective 
potential will spoil inflation, i.e. a slow-roll behavior seems not realistic. We will demonstrate this with a 
region in the {(p, i^)— plane where one still might hope to obtain a sufficiently gentle slope to induce the needed 
long-lasting accelerated expansion of the external space as well as an attraction to the global minimum (in order 
to ensure a late-time stabilization of the scalar fields). (See Figs. 7, 8.) Such a region could be expected close 
to the maximum of the potential U{^) where in rough approximation holds 

U^ff^^^^U {<!>). (121) 

The action functional (20) shows that the two scalar fields ip\ := 0, (^2 := ^ hve in a flat (a— model) target 
space. Hence, the estimate of the slow roll parameters can be performed within a simplified version of the 
multi-field inflation scheme of Refs. [8, 11, 52]. Assuming in rough approximation that the external space has 
already flattened, the inflation parameters e, r] read 

e = -A « (122) 

EU^i^i J:h=MjUeff){diUeff){djUeff) 

Here, H = a/ a is as usual the Hubble parameter of the external space. Additionally, the following abbreviations 
have been introduced 

2 2 

a,:=A, lOUejff = Y,{d,Uefff, \^\' = Y,V>1 (124) 

j=i i=i 

Inflation is possible for e < 1 and a sufficiently small slow-roll parameter \r]\ <^ 1. In the vicinity of the maximum 
in 0— direction it holds dcj)U\max ~ so that e, rj are essentially defined by the slope in (p— direction. Explicitly 
one obtains in this region 

e«2s2, 77«2s(l-s). (125) 
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Figure 8: 3D plot of the effective potential Ueff{4>, (fi) given in Eq.(118). 



For a four-dimensional external space Dq 



and the rough estimates 



rfi = 2 
di = 3 
di=6 



4 this yields from (118) 



2{di + 2) 



1 

<2 



0.5, 


7]^ 


s 0.5, 


0.6, 


77 s 


:^ 0.49, 


0.75, 


7?S 


s 0.47 . 



(126) 



(127) 



Because \rj\ ^ 1 is not satisfied, the considered toy model would produce an accelerated expansion (e < 1) which 
would be much too short for successful inflation. Whether a domain wall with possibility for topological type 
inflation could form between regions to the left and to the right (in 0— direction) of the crest at the maximum 
of C/(^), remains an open question. Corresponding indications have been given in Ref. [19], but seem to require 
an additional detailed analysis — and probably an embedding of the toy model into a more general setup with 
richer structure. 



6 Conclusion 

In the present paper we continued our investigation [33, 34] on multidimensional gravitational models with a 
noii-Einsteinian form of the action. The corresponding action functional was assumed as a smooth function 

f{R) of the scalar curvature _R of a _D— dimensional spacctimc manifold with warped product structure. The 
main subject of our considerations was the stabilization problem for the extra dimensions. As technique we 
used a reduction of the nonlinear gravitational model to a linear one with an additional self-interacting scalar 
field (nonlinearity scalar field 0). The factorized geometry allowed for a dimensional reduction of the considered 
model and a transition to the Einstein frame. As result, we obtained an effective four-dimensional model with 
nonlinearity scalar field and additional minimally coupled scalar fields which describe conformal excitations of 
the scale factors of the internal space (its zero- mode excitations). 

In terms of these scalar fields we performed a detailed stability analysis for models with scalar curvature 
nonlinearities of the type f{B) = R — i-i/R, /it > and /{R) = R + — 2A£), where A/j plays the role 
of a D— dimensional bare (bulk) cosmological constant. As stability condition we assumed the existence of a 
minimum of the effective potential of the dimensionally reduced theory so that a late-time attractor of the 
system could be expected with freezing stabilization of the extra-dimensional scale factors and the nonlinearity 
field. It was shown in Refs. [33, 34], that for purely geomc;trical setups this^is only possible for negative scalar 
curvatures, {R< 0), independently of the concrete form of the function f{R). 
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Four-dimcnsional purely gravitational models with Rr^ eurvature contributions have been proposed recently 
as possible explanation of the observed late-time acceleration (dark energy) of the Universe [14] . In section 4 of 
the present paper, we showed that higher dimensional models with the same scalar curvature nonlinearity 
reproduce (after dimensional reduction) the two solution branches of the four-dimensional models. But due 
to their oversimplified structure these models cannot simultaneously provide a late-time acceleration of the 
external four-dimensional spacetimc and a stabilization of the internal space. A late-time acceleration is only 
possible for one of the solution branches — for that which yields a positive maximum of the potential U {(f)) of 
the nonlinearity field. A stabilization of the internal spaces requires a negative minimum of U ((/>) as it can be 
induced by the other solution branch. The question of whether this incompatibility could be resolved by different 
scalar curvature nonlinearities over the factor spaces (for each factor space Mi it might hold its own curvature 
nonlinearity fi{Ri))\s still open and deserves a separate analysis. We left this issue to future investigations. 

The considered i?''— setup was assumed as highly simplified toy model analogue of the loop corrected gravity 
sector of M— theory [46] . The stability analysis of the higher dimensional model was reduced to a set of algebraic 
compatibility tests for the extremum condition (in the present case a quartic equation in the scalar curvature R) 
and inequalities which ensure the existence of a minimum of the effective potential (non-tachyonic mass terms 
of the corresponding field excitations). For simplicity, we restricted the investigation to parameter regions 
/' = 1 -h 47 > which are smoothly connected with the curvature-linear model at /' = 1 without passing a 
conformal singularity. In Brans-Dicke (Jordan) frame the latter requirement ensures an effective gravitational 
constant which is positive definite and smoothly connected with that of a given BD frame model with a fixed 
(frozen) gravitational constant. 

With the help of a projection technique in the (A/j, 7, i?)— space M we identified regions which ensure 
the existence of a minimum of the effective potential, and hence fulfill a necessary condition for a successful 
freezing stabilization of the extra dimensions. The results can be summarized as follows. For systems with 
total spacetime dimensions D > 6 (in the case of a (Dq = 4)— dimensional external spacetime) there exists a 
stable sector 6i_(a„_^) = {Ad < fl 7 > fl \z{h.D,^)\ < \w{D)\} on the (A^, 7)— plane which in the limit 
7 +0 tends smoothly to the .R— linear sector. The corresponding transition is connected with a freezing of 
the nonlinearity field at the minimum of its potential U{(l)), i.e. a diverging excitation mass, +00, due to 

a diverging Hessian of the potential U {4>). Models within 0i,(Ad.7) are separated from the conformal singularity 
at /' = (and the antigravity sector /' < beyond it) by a potential barrier of infinite height and width and 
are, hence, absolutely stable with regard to transitions into the (/' < 0)— sector. Additionally, it was shown 
that the limit A^ — > —0 in the Gi (A^ ,y)— sector, is connected with a decompactification of the internal space 
components M^, i = 1, . . . , n and a flattening R —0 of the bulk spacetime M. 

Apart from 6i_(Ao,7) there exists a second stability sector 62,(Ar,,7) = {7 < <^ — 1 < z{A£),'y) < —'w{D)} 
for dimensions D < 8. The potential U{(f)) for such configurations is unbounded from below in the limit /'(</> — > 
—00) +0 and has a minimum which is separated from the conformal singularity at /' = (and the antigravity 
sector /' < beyond it) by a potential wall of finite height and width. Configurations in this minimum would 
be metastable and prone to collapse into /' = 0. The 62, (Ad, 7)— sector is disconnected from 9i_(Ao,7) 
an essential singularity of R and U{(j>o) in the limit 7 —0. In this limit the nonlinearity field (j) freezes at 
00 (7 —0) (IM) ln(3r'/|D— 8|) with — > -|-c», but simultaneously the scalar curvature diverges R — > —00 
and the potential deepens unboundedly U{(f)o) —00. This behavior is a strong indication for inconsistencies 
of 02,(A£,,7)~configurations within the framework of the given limited setup. The question of whether the 
62, (Ad, 7)— sector of the considered oversimplified toy-model will find a physically sensible interpretation within 
a still unknown extended curvature-nonlinear theory of gravity, a special UV limit of non-perturbative M-theory 
or within loop quantum gravity remains an open issue. 

A further issue which was out of the scope of the present paper was the analysis of dynamjcal transitions 
between configurations which correspond to different solution branches R{(j)) of Eq. (6), f'{R) = e^'^. For 
the model of section 4, e.g., there exist two such branches R{(j)) which form a double cover over given 

values of the parameter fj, and the nonlinearity field <p. At early evolution stages of the Universe, transitions 
between these two branches cannot be ruled out a priori and should be taken into account for a comprehensive 
description of the dynamics of the Universe. 

Finally, we note that the external spacetimc in the considered pure geometrical models is necessarily AdS 
and the corresponding negative effective cosmological constant, Agjf < 0, forbids a late-time acceleration. The 
situation can be cured by including additional matter fields. Examples are flux field stabilization scenarios [34] 
which provide certain moduli space sectors with positive effective cosmological constant and external spacetimes 

^^See footnote 18 for a few comments on the BD-antigravity sector with /' < 0. 

^^This is in contrast to the considered .R''— model whose solutions R{<t>) are unambiguously defined by the sign of the nonlinearity 
parameter 7. The origin of this difference between the R~^— and the .R*— model is in the number of real- valued solutions Ri{4i) of 
Eq. (6) for given nonlinearity parameters or 7 and 4>. In the case of the model the corresponding quadratic equation has 

two real- valued solutions for given (p) whereas the cubic equation (55) has only one for fixed (7, cp). 
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A The quartic equation and its associated quadratic equation sets 

In this appendix we briefly re-dcrive the quadratic equation sets associated to the quartic equation (59). Fol- 
lowing the Ferrari formalism as it is briefly described, e.g., in [47] we lay explicit emphasis on the sign rules 
which are crucial for a correct derivation of the final solution set of the quartic equation. 
The crux of the Ferrari formalism applied to a quartic equation of type (59), 

+ aix + ao = 0, (128) 
consists in factoring it by transforming it into a difl[erence of two quadratic terms 

A'^ -B'^ = {A + B){A-B) = (129) 

so that solutions can be obtained from 

A±B = 0. (130) 
Adding and subtracting a term x^u + (u/2)^ in (128), with u an auxiliary function, one rewrites (128) as 

(»'n)^-('^-v-^)- 

and requires the second term to be quadratic 



Here, e is a sign factor e = ±1 and we assume for definiteness 



— an „ tU^ — an . s 

i ^ > for i > 0. (133 

u u 



The compatibility of Eqs. (131) and (132) is ensured by the condition 



-^=2eji^^ (134) 
u \ u 

which on its turn is equivalent to the cubic equation 

- Aaou -a\ = Q. (135) 

We specify as in (59), (63) 

^0 ^=~ ^(iP-8) ' ^= 7(J^-8) ' -4^0 = 3«' ^^1 = 2^ (136) 

and assume that u is a real-valued solution of Eq. (135). The analysis of section 5.1 shows that for stable con- 
figurations of the i?^— model it holds additionally u > 0. Using this condition as simplifying input information, 
we can rewrite the quartic equation (128), (132) as 

+ (V^x+|v'w' + 3g)' =0. (137) 



i^These sign rules are not displayed explicitly in [47]. 
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It factorizes according to (129), (130) into a set of quadratic equations 



1 



±^^x + -i^u±e^h?Hiqj =0. (138) 



The sign factor e follows from (133), (134), (136) and u > as 



/ D — 2 

e = — sign(ai) = —sign ( — — — — ) . (139) 



V7(^ 

B Sign analysis of the discriminant Q 

The sign of Q can be obtained by mapping the minimum-ensuring inequality (60) into an equivalent inequality 
for z. For this purpose we consider the critical surface Sci C A4 in the parameter space, where for 7^0 the 
inequality (60) is replaced by an equality, 

E,i = {{AD,7,R)eM\ Cei[AD,7,i?]:=(^-2)+4(D-8)7i?3 = 0}. (140) 

The intersection of this surface Sd with the algebraic variety V : h[Ai:),j,R] = of the extremum condition 
will define a critical value Zd- This value can be found explicitly by resolving (140) for R, what gives 

o \ 1/3 



Rd:=R\_ =-{— — , (141) 

and plugging R^ into the quartic equation (59). As result one obtains 

Zci{^D,l) = 47(8Az3/3)3w(£)) = -1 . (142) 

Now, small perturbations off the critical surface Sci , but along the variety V, can be used to map the inequality 
(60) into its counterpart for z. Setting 

Ad = Ad.ci +<5Ad, R = Rd + SR, z = Zci+5z (143) 

and keeping 7 fixed, we get from the inequality (60) 

72{D - 8)j'^Rl^6R > (144) 

whereas the quartic equation (59) and the definition (66) of z yield 

SAn = S^^^Rli {SRf , Sz = 12^w{D){8/3fAl^dS^D (145) 

and hence 

Sz = (7i?ciAD,ci)^ (SR)' > 0. (146) 

We observe that, although inequality (144) implies 

D>8: 6R>0, D < 8 : 6R < , (147) 

independently of the signs of 7 and -Rci, the variety V is for 7 ^ 0, ^ 8 and (because of 6z > 0) located over 
the region 

^(Ad,7)>^ci(Ai3,7) = -1 (148) 

of the (Ad, 7)— plane. This means that by any perturbation (motion) on the variety V we cannot pass across 
the critical value Zci(Ad7 7) = — 1- Hence, Zci(Ad,7) = — 1 must be a boundary segment of the projection ttV 
of V onto the (A^, 7)— plane: z^i C 9(7rV). The latter fact is confirmed by the observation that the critical 
surface Sd C A4 coincides with the singular surface 

a^/i[AD, 7, R] = {D-2)+ 4{D - 8)7i?3 = 

of the projection tt of V onto the (Ad, 7)— plane. 
^''Singularities of smooth projections are extensively discussed, e.g., in Refs. [53]. 
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C Mapping f > into parameter space 



The inequality /' = 1 + i'yR^ > can be analyzed with the same technique as the minimum-ensuring inequality 
(60) (see relations (140) - (148)): we obtain the intersection of the critical surface 

Ec2 = {{Ad,J,R)&M\ U[^d,1,R] = 1 + ^jR' = 0} (149) 

with the algebraic variety V, i.e. Sc2 n V, and study the behavior of small parameter perturbations off Sc2 and 

along the variety V. 

Explicitly this means that we resolve (149) for R to obtain 

Rc2 = -(47)"'/' (150) 
and plug this Rc2 into the quartic equation (59). From the intermediate result 

47(8Ab/3)3 = -1 (151) 
we find by multiplication with 'w{D) that the intersection Sc2 H V corresponds to the critical value 

Zc2{Ad,i) = -w{D). (152) 
Substituting, furthermore, the perturbation ansatz 

An = Ad,c2 + SAd, R^Rc2 + SR, 
z = Zc2 + Sz, f' = Sf, (153) 

(it holds = 0) into the defining relation (66) for z, the quartic equation (59), and the equality /' = 1 + 47^^, 
we get for fixed 7 

5z = l2-fw{D){%/^fK^jj^^^5KD, 6Ad = ^6R, 

6f = UjRl^SR, (154) 
respectively, and by combination of these results also 

1 ^ ./ ,- 3 /8\ / Ai3_c2 ^ ^ 



Prom the definition (66) of w{D) and its implication 

w{D < 8) < 0, < w{D > 8) < 1 (156) 

we find for 6f'>0, f > 

D<8: Zc2{Ad,i) = \w{D)\ > 0, z< -w{D) = \w{D)\ , (157) 

Z)>8: Zc2{Ad,i) = -w{D) <Q, -w{D) < z . (158) 



^*For completeness we note that a formal crossing of the critical surface Sc2 into the sector /' < would lead to the restrictions 

D<8: z>\wiD)\, D>8: z < -w{D). 

The sector /' < itself would correspond to a negative effective gravitational constant in BD frame (an antigravity sector) [54, 55]. 
Antigravity effects arc described for various types of higher dimensional setups. In certain SUGRA configurations, spatially bounded 
antigravity regions arc known as repulsons (white holes) [56] . Several brane world models show radion induced scalar antigravity at 
ultra-large distances [57]. In our model, the antigravity sector would fill a complete external spacetime which necessarily would be 
disconnected from our own observable Universe with /' > by a conformal singularity at /' = (see subsection 5.3 below and also 
[55] ) . There remains the interesting open question of whether the formal correspondence between the strong-curvature sector of the 
^—nonlinear model and the emerging antigravity sector of the R— linear model in BD-frame simply signals an inconsistency of the 
theory and whether it can find a physically meaningful interpretation in an enlarged (extended) setup. The repulson in SUGRA 
models was reinterpreted as unphysical region and found its resolution by the enhangon mechanisms (excision of the antigravity 
region and placing a heavy shell of wrapped D-branes on its boundary so that the former antigravity region in the interior is shielded 
by the D-brane source and replaced by a segment of flat spa<;e) [58] . In analogy to the enhangon mechanism in the repulson case, 
one can expect a resolution of the conformal singularity at /' = by some kind of quantum gravity mechanism. 
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D Parameter limits 



In subsection 5.3 it has been found that for the (/' > 0)— sector the boundary segments ^(A^)."/) = —w{D) C 
90(Ac,7) of the projection 6(Ad,7) ■~ '^'^ of the stabiUty region T cV C A4 onto the (A£),7)— plane correspond 
to the limit (j) — * — oo. Here, we clarify the behavior of the system in the vicinity of the other boundary segments 
dO(Ann) ^ i^D = U 7 = 0} for £) > 8 (see Fig. 4) and Se(Ac,7) ^ i^D = OU^ = OU z = -1} for Z? < 8 
(see Fig. 3). 

Ad ^ 0,7 ^ : In this limit we obtain from Eqs. (63), (66), (68) 

q^O, z^O, Q ^ ^ 0, vi ^ 2^/^ (159) 
and hence from Eqs. (76), (77), (92) 

7 > : ^(Ad ^ -0) ^ -0 ^ { R ^ -0 for D = 8, (160) 



R-.+ ^ 


-0 


for 


D > 


R 


-0 


for 


D = 




-0 


for 


D < 


for £> < 8 









7<0: i?_ _(Au ^ 0) ^ -(2r)^/*^ for £> < 8. (161) 

Obviously, the system behaves differently in the upper and lower (A^j, 7)— plane. In the case of 7 < 0, the 
system behaves regularly for A/) and the half- line (A^) = 0,7 < 0) is not distinguished from its vicinity. 
In contrast to this, the limit (Ad ^ —0,7 > 0) corresponds to a flat-space limit ii — > — which via (55) is 
associated with a freezing of the nonlinearity field (j): /' — > 1 at (po — > 0. From Eq. (41) follows 

so that for the mass of the nonlinearity field holds to| — ^ +00. We note that the nonlinearity field in an 
i?^— model has a finite mass in the limit A^ —0 (see, e.g., [33, 34]). The different behavior of the models 
is caused by the different powers of the term {e^"^ — 1) in U{(j)): for an i?^— model this power equals 2, whereas 
for an i?^— model it equals 4/3. Hence, in the latter case the second derivative (PU{(j))/d(j)'^\^^^^ diverges in the 

limit (/)(o) ~^ 0. 

We arrived at the interesting fact that in the considered toy model the extremum condition in form of the 
quartic equation (59) relates the scalar curvature R at the minimum and the bare cosmological constant A^ in 
the case of 7 > so strongly that for stabilized internal spaces the limit A^ — > corresponds to the fiat-space 
limit R — > —0. As it should be, the flat-space limit of the total scalar curvature R —0 implies via (43), i.e. 
U{(j)o) —0, and (32), Ri = 2diU {(j)o) / {D — 2), also a decompactiflcation of the internal space components 
Ri = e-'^^'oRi -0, +00 (the Ri are hold flxed). 

7 — »■ 0, Ad 7^ : The deflnitions (63) and (66) show that for fixed A^ ^ the limit 7 — > implies 

r,\q\,Q^ +00, z^O. (163) 

With the help of an expansion in terms of small z ~ the rescaled curvatures (77) are easily obtained from 
(69) as 

r+_(z^0) Ri -3 X 2-^/2^;i/3 
T_,+(^^0) « 3X2-5/V/3 

T_,_(^^0) w _2V6 _ 2-5/2^1/3 (164) 
so that the curvatures Re,± themselves can be estimated via Re,± = r^^^Te,± as 

^_,+(7^+0;£>>8) \ _ 2£>Aj, 



i?+,_(7^+0;£><8) j £>-2' 



i?_,_(7^-0;£)<8) 



D-2 



^{D - 8) 



(165) 

1/3 

^-00. (166) 



^^The general limiting behavior R(7 — > 0) without identification of the concrete solution branch fle,± can be easily obtained from 

the quartic equation (59). Assuming _R(7 — » 0) < 00 and taking the limit 7 — > in Eq. (59) gives R = 2DKo/{D — 2), whereas 

1/3 



division of (59) by R for a behavior |-R(7 — » 0)| — > 00 yields R = — ^ ^l^_^g'^ ^ 
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In the exceptional {D = 8) case the scalar curvature in the minimum does not depend on 7 and is given by Eq. 
(92) 

i?=«A, = ^. (167) 
3 D-2 ^ ' 

Again, the system behaves qualitatively different in the upper and the lower (A£), 7)— plane. Because of the 
finite asymptotics (165) and Eq. (167), in the upper half-plane it holds (for D < 8 and D > 8) 

/'(7 ^ +0) ^1, (l>o^0 (168) 



f/('^o;7^+o) 



{D - 2)3 



96{D- 
Ad 



+00 



(169) 
(170) 



and the nonlinearity field (/) undergoes a freezing stabilization at 00 = with diverging mass +00 but 

finite scalar curvature (165) and finite minimum position U{4)o). Hence, under the freezing stabilization of the 
nonlinearity field for 7 +0 the system turns srnoothly into a system with linear scalar curvature term R, 
i.e. into a system with Einstein-Hilbert action in R. This is a generic feature of models with nonlinear scalar 
curvature terms and was earlier described for i?^— models in [33, 34]. Figure 5 gives a rough illustration of the 
corresponding deformation of the potential U{(j)) under variation of 7 and for a fixed value of A/j. 

The behavior of the system is completely different in the lower- half-plane limit 7 — > — 0. Here we have to 
distinguish the dimensions D = 8 and D < 8. For D = 8 Eqs. (169), (170) extend to the lower half-plane 
7 — ^ —0, i.e. the system is completely unstable in this limit d'^U\ct,„(j —0) —00. Obviously, d'^U\4,o ^ 1/j 

in (169) encounters a pole singularity with respect to 7. This is different for R , D < 8. Here, one finds from 

(55), (166) 

3D 



= /'(7 - -0) 



and from (41), (43) 



1 (Z?-2)i/3(i?- 8)2/3 



^1/3 
1 {D 



8{D 

2)4/3 



71/3 2D\D-8\) 



-1) 
3D 



3D 



D 



-|-oo 



D-8 



—00. 



(171) 

(172) 
(173) 



Prom these equations and their rough illustration in Fig. 6 we read off that in the limit 7 — > — the minimum 
of the potential ?7(0) lowers infinitely, U{(f)o;j — > —0) — > —00, and becomes fixed at the finite value 0o = 



3D 



D-8 



The corresponding scalar curvature diverges as i?_,_(7 —0) ^ —00 and is separated from the 

conformal singularity /'(0 — * —00) ^ by a barrier whose top is defined by the associated maximum branch 
R-,+ and tends to the value 



^_,+ (7^-0) 



2DKd 
D-2' 



(7 ^ -0) ^ Kd 



(174) 



(see also Fig. 6). This means that the metastable sector in the lower (A/?, 7)— plane is separated from abso- 
lutely stable systems in the upper half-plane and their limiting linear (7 -|-0)— models by an infinite gap. 
Considering the behavior of the system over the (A^), 7)— plane we have to conclude that it possess an essential 
singularity in the limit 7 —0, i.e._ the scalar curvature R encounters an infinite jump between 7 +0 
[R 2DAd/{D — 2)) and 7 ^ — (i? — > —00) whereas 9^f7|<^o for 7 > a pole-like singularity ~ 7"^ and 
for 7 < it behaves like a branching point singularity ^ 7"^/^. This is a strong indication that the description 
of a physical system in terms of the considered oversimplified toy model breaks down in the limit 7 —0. The 
search for possibly existing physically realistic metastable systems with 7 < is out of the scope of the present 
work and we leave a corresponding investigation to future research. 

z —>■ —1 :_Relations (88) - (91) show that this limit can only be reached by metastable configurations 7 < 
0, D < 8, i?_ _. According to (66) it corresponds to a vanishing discriminant Q = of the cubic equation (62) 
so that this equation has two coinciding solutions Mi, 2 = 2r^/3 holds vi^2{z = —1) = 2). From Eqs. (76), 
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(77) we observe that this leads to a coalescence of the minimum branch i?-,- and the associated maximum 
branch -R_,+ in an inflection point at 2; = — 1 

= -1) = -ly^J^ = -r!^ = -2-»/= (^) . ,175, 

The situation is also obvious from Fig. 2. Configurations with z < —1 have no extremum at all and are 
necessarily unstable. 
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